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ABSTRACT 


The  long-range  forces  between  large  conjugated  double  bond  molecules  are 
highly  direct  ional  and  anomalously  strong.  This  is  due  to  pi -electrons  moving  free- 
ly from  one  end  of  the  conjugated  double  bond  system  to  the  other,  making  the  mole- 
cule behave  tike  an  emended  oscillator  All  three  types  of  dispersion  forces:  CT  cr 
CTT  and  7T  TT  are  calculated  for  ethylene,  acetylene,  and  the  higher  polyenes.  In 
addition  the  Tr  TT  interaction  energy  is  calculated  for  benzene.  The  free  electron 
molecular  orbitals  are  used  for  the  TT  electron.^  and  products  of  these  orbitals 
are  approximated  by  "transition  monopoles"  Our  treatment  of  the^^electrons  and  the 
use  of  the  transition  monopoles  were  suggested  by  London^  The  resulting  energies 
are  easy  to  calculate  and  easy  to  interpret.  The  TfJT  energy  is  proportional  to  the 
square  of  the  energy  of  interaction  of  two  real  transition  dipoles  (real  in  the  sense 
of  the  monopoles  being  separated  by  a finite  distance).  Previously,  Coulson  and 
Davies^  calculated  the  interaction  energy  using  LCAO  orbitals  and  evaluating 

all  integrals  accurately,  however,  because  of  the  complicated  nature  of  their  cal- 
culations it  was  not  possible  for  them  to  obtain  a simple  understanding  of  the  varia- 
tion of  the  energy  with  the  orientations  of  the  molecules.  The  values  we  obtain  for 
the  '/ITT  energy  agree  with  those  of  Coulson  and  Davies  except  for  a scale  factor, 
our  results  being  uniformly  3/5  .small.' r.  Such  a .scale  factor  corresponds  to  an 
adjustable  parameter  in  either  the  Coulson  and  Davies  treatment  or  ours.  The 
CrTT  interaction  energy  has  never  previously  been  considered.  For  long  polyenes, 
for  which  a simplified  treatment  is  given,  the  V TT  energy  is  larger  than  the  O'  F 
which  is  larger  than  that  of  the.  (7(7  . In.  the  appendix, the  free  electron  model  is 

used  to  calculate  the  pi-electron  contribution  to  the  polarizability  and  excellent 
agreement  is  obtained  with  the  LCAO  calculations  of  Davies^. 
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I.  Introduction. 

The  long-range  dispersion  forces  between  large  conjugated  double  bond 
molecules  are  highly  directional  and  very  strong  This  is  due  to  pi-electrons 
moving  freely  from  one  end  of  the  conjugate  double  bond  system  to  the  other, 
making  the  molecule  behave  like  an  extended  oscillator.  These  forces  may  help 
to  explain  some  of  the  anomalous  effects  of  steroids  and  other  conjugated  double 
bond  molecules  of  biological  importance,  Also  the  pi-bond  forces  are  important 
in  the  interaction  between  aromatic  molecules  and  are  partially  responsible  for 
the  non - zero  angles  between  the  planes  of  neighboring  molecules  in  an  aronaatic 
crystal.  In  this  paper,  a simple  treatment  is  developed  for  calculating  and  ex- 
plaining these  forces. 

In  contrast  to  the  pi -e  lectrons,  molecules  contain  sigma-electrons 
which  are  localized  in  specific  chemical  bonds  or  within  the  inner  shells  of  the 
atoms.  Corresponding  to  the  two  types  of  electrons,  the  dispersion  energy  can 
be  divided  into  three  terms; 

-H  n.l) 

Here  E (j- q-  is  that  part  which  arises  from  the  interaction  between  the 
(T  - electrons  of  the  two  molecules;  E corresponds  to  the  interaction 
between  the e lectrons  of  the  two  molecules;  and  E^— 7^  gives  the  energy  of 
interaction  of  the  — electrons  of  one  molecule  with  the  7K—  electrons  of  the 
other  molecule  and  vice  versa.  London^ developed  an  excellent  method  for  cal- 
culating the  Eq-(j-  • It  is  the  purpose  of  the  present  paper  to  describe  a simple 
method  for  calculating  E (j- 7^  and  E -p- 7^  • 

We  use  free  electron  molecular  orbitals  for  the  pi -electrons  and  approx- 
imate products  of  these  orbitals  by  "transition  monopoles".  The  resulting 


1.  E'.  London,  J.  Chem.  Phys  , 46,  305  (194<i| 
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values  of  E q— ;7-  and  E are  easy  to  calculate  and  easy  to  interpret. 

For  example,  to  a good  approximation  in  most  cases,  the  values  of  Ep-;^ 

are  proportional  to  the  square  of  the  energy  of  interaction  of  two  real  trans  • 

ition  dipoles  (real  in  the  sense  that  the  monopoles  are  separated  by  a finite 

distance).  We  are  indebted  to  Eondon^  for  the  notion  of  transition  monopoles 

and  their  use  in  connection  with  dispersion  energy.  Our  results  for  E^^ 

2 

are  to  be  compared  with  the  calculations  of  Coulson  and  Davies  who  used 

LCAO  orbitals  and  evaluated  all  of  the  integrals  accurately  without  using  the 

monopole  approximation.  Because  of  the  complications  of  their  calculations 

it  was  not  possible  for  them  to  obtain  a simple  understanding  of  the  variation 

of  the  interaction  energy  with  the  orientations  of  the  molecules.  However, 

our  values  of  agree  with  those  of  Coulson  and  Davies  except  for  a scale 

factor  (our  values  being  uniformly  3/5  of  the  Coulson  and  Davies  numbers). 

Such  a scale  factor  corresponds  to  an  adjustable  parameter  in  either  treatment 

Prior  to  the  present  paper,  there  were  no  evaluations  of  E 0-77-  . For  long 

polyenes  by  no  means  negligible  and  as  a matter  of  fact:  . 

For  ethylene,  acetylene,  and  the  lower  polyenes:  E^^^y  E.^^')Eyf.^  . In 

the  case  of  ethylene  and  acetylene  where  the  complete  dispersion  energy  is 

compared  with  experimental  data,  the  agreement  is  satisfactory  in  the  sense 

that  it  is  quite  comparable  with  the  ratio  of  theoretical  to  observed  values  for 

simple  molecules  in  which  there  are  no  complications  due  to  pi— electrons. 

In  the  appendix,  the  free  electron  molecular  orbitals  are  used  to 

calculate  the  contributions  of  the  pi-electrons  to  the  polarizability  of  polyenes 

3 

Excellent  agreement  is  obtained  with  the  LCAO  calculations  of  Davies  and 
the  results  seem  to  be  consistent  with  experimental  data.  This  agreement 
serves  as  a check  on  the  accuracy  of  the  free  electron  molecular  orbitals 
as  used  in  the  present  paper. 


2.  C.  A.  Coulson  and  P.  L.  Davies,  Trans.  Faraday  Soc.  48,  777  (1952), 
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In  retrospect  we  feel  that  the  simple  free  electron  molecular 

orbitals,  as  used  here,  give  roughly  the  right  values  for  and 

However,  the  fart  that  the  orbitals  have  been  taken  as  one -dimensional 

rather  than  as  three-dimensional  has  its  limitations.  For  example,  consider 

the  interaction  of  two  ethylene  molecules  in  which  the  C = C axes  of  the  two 

molecules  are  parallel.  With  the  present  free  elect  ron  orbitals,  the  value 

of  IS  independent  of  the  orientation  of  the  planes  containing  the 

hydrogen  atoms.  Such  a result  is  obviously  erroneous  An  obvious  im  - 

provement  would  result  by  replacing  our  present  one -dimensional  molecular 

orbital;®  (x)  , 'see  Eq-  2.29)  by  three  -dimensional  free  electron  molecular 
' n 

orbitals  such  as 

^ y (1.2) 

Here  x is  the  cartesian  coordinate  along  the  axis  of  the  conjugated  double 
bond  system  and  y and  7.  ai  e coordinates  perpendicular  to  this  axis.  The 
parameter  oC  , would  of  course  be  adjustable.  The  use  of  three  - dimen- 
sional free -electron  orbitals  would  not  appreciably  increase  the  amount  of 
work  required  and  at  the  same  time  it  should  improve  the  comparison  of 
the  dispersion  energy  for  different  orientations  of  the  molecules. 


2.  1 

11.  Mstbod. 

la  dealing  with  dispersion  forces,  the  molecules  axe  sufficie'otly  separated 
to  a pertiirbation  calculation  correct.  Furthermore,  it  is  not  necessary  to 

antisymmstrixe  the  wave  function 'fof  the  two  molecule  system,  since  we  may  as- 
sume that  there  is  no  overlapping  between  the  wave  functions  of  different  mole- 
cules. We  may  then  employ  a simple  product  of  the  electronic  wave  functions  of 
the  individual  molecules  as  our  zero  order  wave  function  for  the  system  and 
treat  the  potential  energy  of  interaction  as  the  perturbation.  Actually,  the  only 
part  of  the  potential  energy  of  interaction  which  gffects  the  dispersion  energy  are 
the  electron-electron  interaction  terms  '■ 


t-tt 

^ »'<  L--1 


e , 

Tab 


(2.1) 


Here  the  n^  electrons  designated  by  the  subscript  "a"  are  in  molecule  "A" 

and  the  n^  electrons  designated  by  the  subscript  “b"  are  in  molecule 

Furthermore  r . is  the  instantaneous  separation  between  electrons  *'a-' 

A 

”b”.  If  Sto  and  aure  the  electronic  wave  functions  of  molecule  A 

B B 

in  its  ground  state  and  l-th  excited  state  respectively;  if  and 

are  the  wave  functions  of  molecule  B in  its  ground  state  and  j-th  excited 
state  respectively;  and  if  , E^  , and  E^  are  the  energies 

of  these  xnolccular  states  respectively,  the  usual  second  order  perturbatioii 
theory  gives  the  energy  of  dispersion; 


E 


(dis) 


(Z.  f) 


1.  The  matrix  components  which  give  the  greatest  eon.tribntion  to  the  dia° 
persion  energy  correspond  to  one  electron  in  each  of  the  two  noolecules 
making  a transition.  For  such  two  electron  transitions,  only  the  electron - 
electron  terms  in  the  perturbation  potential  as  given  by  Eq.  (2.  1)  lead  to 
non-vanishing  matrix  components.  However,  it  should  be  pointed  out  that 
one  elsctxon  transitions  do  give  small  but  finite  contributions  to  the  disper- 
sion energy.  For  one  electron  transitions*  it  is  the  electron-other-nuclcus 
terms  in  the  perturbstion  potential  which  are  important.  In  the  present 
paper  we  ignore  the  small  contributions  from  one  electron  transitions'.x 
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H»re  the  prime  over  the  summation  sign  indicates  that  those  states  with  either 
or  -fe  equal  zero  are  excluded.  The  conjugate  complex  of  the  wave  function 
does  not  appear  in  Eq.  (2.2)  since  the  wave  functions  in  the  following  d evelop- 
ment  are  taken  to  be  real.  Furthermore,  for  present  ptxrposes  it  is  convenient 
to  suppose  that  the  wave  function  for  each  of  the  molecules  can  be  taken  as  a 
Slater  determinant  made  up  of  molecular  orbitals,  ^ and  0 The 

orbitals  are  numbered  by  a subscript  according  to  their  energy,  6 . Thus 

(p  ^ is  the  orbital  of  lowest  energy,  € ^ ^ 2 next  higher 

energy  , c _ ; etc.  the  molecules  which  we  shall  consider  have  an  even 
number  of  electrons  n = 2N,  and  their  groxmd  state  is  a singlet.  To  the  ap- 
proximation of  the  present  treatment,  there  are  no  singlet-triplet  transitions 
80  that  the  only  excited  states  which  need  be  considered  are  themselves  sing- 
lets. Thus  the  wave  function  for  the  ground  state  of  a znolectile  is  a Slater  de- 
terminant with  two  electrons  in  each  of  the  molecular  orbitals:  0 , • <t  ^ 

wave  functions  of  the  excited  states  which 
we  consider  asze  the  same  as  the  wave  function  of  the  ground  state  except  that 
one  of  the  two  electrons  in  the  i-th  molecular  orbital  has  been  promoted  to  the 
j-th  molecular  orbital.  Such  a transition  is  denoted  by  the  symbol  (i«  j). 

In  terms  of  these  orbitals  the  ground  state  of  the  molecule  A is  the  Slater  de- 
terminant. 


<j>’' ■ ■ ■ (4.''  <6''  gs"  izS"  jzi''  \ 


/9  «.  y5 


) 


(2.  3) 


Here  the  i-th  orbxtals  are  moved  from  their  normal  positions  in  the  deter- 
minant in  order  to  facilitate  the  algebraic  manipulations  which  follow.  Simi- 
larly the  excited  state  A of  the  molecule  corresponding  to  the  transition 
fi«  j)  has  as  its  wave  function  the  difference  of  two  Slater  determinants. 


fZ.4) 
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Now  integrating  the  product  of  ^ and 

coordinate  and  spin  apace  for  electrons  1,2,  . . 
leads  to  the  result. 


sp;  over  all  of  the 
(a-l),  (a+1) 2N^ 


J 'f'o' 


— — (bi  (o.)(h^^(a)(tK(<iU(o.) -hfiCa-lfiCo.) 

VT  ^ ^ 


(2.5) 


Since  all  of  the  electrons  are  equivalent,  substitution  of  Eqs.  (2.  1)  and 
(2.5)  into  Eq.  (2.2)  leads  to  the  following  expression  for  the  dispersion 
energy  in  terms  of  the  molecular  orbitals. 


(2.6) 


Here 


(2.7) 


Sphere  the  integrations  are  taken  over  the  coordinate  space  of  electrons  a 
and  b. 
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In  a molecnle  containing  conjugate  double  bonds  there  are  two  kinds  of 

molecular  orbitals . First,  there  are  the  pi-orbitals  in  which  the  electron 

is  not  localized  in  either  a particular  bond  or  on  a particular  atom  but  rather, 

the  pi-electron  is  free  to  roam  over  the  length  of  the  conjugate  double  bond 

system.  The  other  electrons  in  the  molectiles,  which  we  shall  loosely  refer 

to  as  - - electrons,  are  localized  in  either  individual  atoms  or  chemical  bonds. 

The  usual  treatment  of  dispersion  energy  applies  to  interactions  between  the 

sigma  type  electrons  in  one  molecule  and  the  sigma  electrons  of  the  second 

2 

molecule,  E^o-.  In  calculating  E ^ , the  terms  e are  expanded  in 

inverse  powers  of  R , the  separation  between  the  molecules.  This  expansion 
converges  rapidly  provided  that  the  linear  dimensions  of  the  electron  orbitals 
are  small  compared  to  the  separation  b etween  the  molecules.  For  pi -electrons, 
the  size  of  the  orbitals  are  comparable  with  the  size  of  the  molecules,  so  that 
this  type  of  expansion  would  not  converge  rapidly  when  the  molecules  are  separ- 
ated by  one  or  two  molecular  diameters.  Thus  a somewhat  different  treatment 
is  required  in  the  calculation  of  the  interaction  of  the  pi-electrons  of  one  mole  - 
cule with  either  the  pi-electrons  of  the  other  molecule,  E -p.  ^ , or  the  sigma 
electrons  of  the  other  molecule,  E ^ ,t  . 


The  London  Approximation 

2 

F.  London  suggested  a method  for  approximating  the  matrix  elements  of 

Eq.  (2  7 ) corresponding  to  either  pi-pi  or  sigma-pi  electron  interactions  in 

2 

which  cases  the  usual  treatment  (expansion  of  e in  inverse  powers  of  R) 


2.  F.  London,  J.  Chem.  Phys.  305  (1942) 
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does  not  apply  for  the  values  of  the  intermolecular  separation  which  have  the 
greatest  physical  interest.  First,  he  considered  the  matrix,  element  to  be 
electrostatic  energy  of  interaction  between  two  electrical  charge  distributions: 

A j A 

— (/^  0 <P  J r«.) 

(2.8) 


f.y  (b)  = e 4>P(^) (i>) 


The  matrix  component  can  then  be  written; 


(2  9) 


These  charge  distributions  are  negative  in  some  regions  of  space  and  positive 
in  others,  depending  upon  the  signs  of  the  molecular  orbitals.  Thus  the  nodal 
surfaces  of  ji  (ij)  divides  all  space  into  a number  of  cells  or  zones  k. 
Similarly  the  nodal  surfaces  of  ^ (i’j*)  divides  space  into  a set  of  cells  or 

zones  labelled  k'.  The  noatrix  component  can  then  be  written  in  the  form; 


(2.  10) 


Throughout  each  cell,  the  charge  distribution  maintains  either  a positive  or 
negative  sign.  London  now  defines  a set  of  monopoles  corresponding  to  the 
charge  distributions.  The  charges  on  the  monopoles  are 


(2.11) 


The  monopoles  are  considered  to  be  located  at  the  points  R (k)  and 
B 

given  by  the  first  moments  of  the  charge  distributions  within  their 
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respective  zpne. 


(2.12) 


The  London  approximation  then  consists  of  replacing  the  individual  integrals 
of  Eq.  (2.1  ) by  the  electrostatic  interaction  of  the  corresponding  monopoles. 


(2. 13) 


Since  the  London  approximation  forms  the  basis  for  the  present  calcula* 
tions,  it  is  import  ant  to  examine  its  validity.  Let  us  consider  one  of  the  in- 
tegrals of  Eq.  (2J0) 

Take  R..(k)  to  be  the  center  of  the  coordinate  system  for  electron  *'a*’ 

B 

and  R (k')  be  the  center  of  the  coordinate  system  for  electron  "b*'.  Let  the 

^ ^ /' 

Z-axis  for  both  coordinate  systezns  point  in  the  line  from  R..(k)  towards 

B _ij 

R.,.,(k').  Then  are  the  spherical  coordinates  of  electron  *'a''  and 

: j • 

(r’,6^^')  arc  the  spherical  coordinates  of  electron  ”b"  in  the  two  separate 

3 

reference  frames.  Then  we  can  expand  a bipolar  expansion: 


I ^ ^ \m-  IrnI 

7 -IT  X Z (cose)?  (cose') 

at  n-  0 n'--o  > 


cm  ( (j)  ') 


(2.  14) 


Here  the 


m 


n,  n' 


are  the  bipolar  expansion  coefficients, 


and  “ j £ 


A 

‘j 


3.  See  ‘'Molecular  Theory  of  Gases  and  Liquids*'  by  J.  O.  Hirschfelder, 

C.F. Curtiss,  and  R.B. Bird  (John  Wiley,  1954),  p.834,  Eq.  (12.1-28). 
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If  the  k-th  zone  of  molecule  A and  the  k'-th  zone  of  molecule  B do  not 
overlap  or  if  r = r'  <(}^,  the  bipolar  expansion  coefficients  have  the  very 
simple  form. 


In 


n+rn  ri  V' 

R ■ / ' r- r r 


(2.  15) 


In  this  case,  the  integral  of  Eq.  (2.10)  becomes: 


f 

J J 
^k)(k') 


Aty? 

1 ® i ' 


JH  -'j 


v/here  the  K ^ are  a set  of  constants  given  by 


s s_  s; 


n '"t-n 


^Aor(s-iyz 


l-rri 


(2.16) 


(2.  17 ) 


The  summation  going  to  J./2.  or  ( I- I )/z  means  to  take  whichever  of  the 
two  is  an  integer.  The  constants  in  this  equation  are  the  integrals^ 

\ ( n . 

. *"  1 ' r . 

n nl  J n 


.1  \ r 1^^(cos6)  § Ck)cl^; 

••J 


a. 


(k) 


s 


(-1) 


i-r)  Q-n)l  J, 
(k) 


r / 

\ p ( cob  o')  § ,(k ) dt 

J ^ i-V  ^ 


(2. 18) 


(2.19) 


L 


r,  (ntm 


)!  ) ^ 0)CO5  >A^  9..(k)c)^ 


(k) 


(2.20) 
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L 


(-1)  ( ;(A-n) 


i-n  (i-n  + no)l. 


6 . 

) '('  (cr.5  ejCOSrm^'?,  / k ) d 


(2.21) 


w 


n,nn  (n+m 


n T A 

( CO^  6 ) sin  no  c;^  §,.  (k)  d 


(k) 


(2.22) 


w 


(-1) 


i-n,rY)  (j[-^  + iY>')|  J 


(k') 


P (cos  0 )Sir)  ? (k)  dt 

(1-n)  t 


A B 

Thus  S - 0..  (h),  S'  =C  ...,(k’)  and  because  of  the  location  of  tlie  origins  of 
o c'lJ  o 'I'l'  'W„-w;, 

the  two  coordinate  systems,  S,  = S'  = L,  , = L'  , = 0.  As  a result  the  lead- 
• 1 il 

ing  terms  of  Eq.  (2.  13)  are 


tj  b 


(2.24) 


0^ . 


4? 


The  London  approximation  therefore  consists  of  taking  the  first  term  of  this 
infinite  series.  The  second  term  in  the  series  vanishes  because  of  the  loca- 
tion of  the  monopoles.  The  third  term,  which  gives  the  first  corrections  to 
the  London  approximation  depends  on  the  second  moment  of  the  charge  distri- 
butions. The  series  converges  rapidly  and  the  London  approximation  is  accur- 
ate provided  that  X is  large  compared  to  the  stun  of  the  effective  radii  of 
the  two  charge  distributions  in  the  particular  zones  being  considered. 


The  Coulson-Davies  Method. 
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Instead  of  using  the  London  approximation,  Coulson  and  Davies  expressed 
the  molecular  orbitals  for  pi-electrons  as  linear  combinations  of  the  2p  atomic 
orbitals, 


(a) 


(Z.25) 


Here  r is  the  distance  of  the  "a"-th  electron  from  the  oc-th  atom 

(carbon  atoms  in  the  problems  considered  here);  the  J?-axis  is  taken  per- 

pendicular to  the  axis  of  the  molecule  so  that  is  the  coordinate  of  the  a-th 
electron  in  this  direction;  and  "c”  is  a constant  which  for  carbon  has  the 
value  1.625/a^  where  is  the  first  Bohr  radius.  In  a linear  conjugated 
polyene  with  2N  carbon  atoms,  there  are  2N  molecular  orbitals  that  can 
be  formed  by  linear  combinations  of  the  above  atomic  orbitals; 


3.N 


od“  I 


(a) 


r 


2N 


(2.26) 


The  one  electron  energies  associated  with  these  orbitals  are 

j = 1,  2,  ....  2N  (2.2Z) 

Here  Y is  the  "resonance  integral"  which  C oulson  estimated  empirically 
to  have  the  value  -40  kcal/mole. 


= 2y  cos  (^) 


4.  C.  A.  Coulson  (Proc.  Roy  .Soc.  169A,  413  (1939))  developed  the  molecular 

orbitals  given  in  Eqs.  (2.22)  and  (2.23)  and  C.  A. Coulson  and  P.L. Davies 
(Tr.  Faraday  Soc.  ^ , 777  (1952))  applied  them  to  problems  involving  pi-bond 
intermolecular  forces. 
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Using  the  moleculaT  orbitals  of  Eq.  (2.  26)  the  matrix  elements  of  Eq.  (2.7) 
which  occur  in  the  dispersion  energy  become  linear  combinations  of  integrals 
of  the  type. 


A ■-  'j  - . ^ ’ 

fi  ^ ^ J -'. 


(2.2&) 


Because  of  the  exponential  decrease  of  the  atomic  orbitals  with  distance,  the 
only  integrals  which  have  .ippreciable  size  correspond  to  c<  and 
cither  coinciding  with  or  being  next  door  neighbors  to  evi  and  /j  respectively, 
Furthermore,  since  the  size  of  the  atomic  ortibals  is  small  compared  to  the 
separation,  R,  between  the  molecules,  it  is  possible  to  evaluate  integrals  of 
the  type  of  Eq.  (2.  ) in  the  same  manner  as  is  usually  done  for  the  matrix 
elements  of  Eq.  (2.7)  in  the  case  of  sigma-sigma  electron  interactions. 


The  FEM  orbitals 

Recently  a different  type  of  molecular  orbital  has  been  introduced  which 
is  more  satisfactory  for  the  calculations  that  follow  because  of  its  greater  sim- 
plicity.  Moreover,  it  is  in  no  way  inferior  for  predicting  the  energies  of  ex- 
cited states  and  contains  no  constants  which  require  empirical  determination. 

This  approach  is  known  as  the  free  electron  model  approximation.  The 
delocalization  of  the  pi-c  le  e irons  is  taken  literally  and  they  are  regarded  as 
free  to  move  in  a one  dimensional  box  which  extends  along  the  skeleton  of  the 
conjugated  carbon  atoms.  Thus  the  skeleton  for  a benzene  molecule  is  essen- 
tially a circle  around  the  benzene  ring  and  for  a linear  polyene  it  is  a line 
along  the  carbon  netwo-k.  Whenever  there  is  a free  endpoint  such  as  the 
terminal  carbon  atom  of  a linear  polyene,  the  ’’box*’  is  extended  an  additional 
bond  length,  or  the  ca.rbon-rarbon  distance. 

For  linear  polyenes  with  N conjugated  double  bonds,  the  free  electron 
model  molecular  orbitals  are  given  by: 


. fx  . niiX.  , ^ 
0 ”'\/T  Siri n / ,2  , 


. . . '.  ♦ / 
u \ y 


(2.2S) 


-I- 
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y* 


with  the  energies: 

r * ^ 

= — Ti  , o- 1,  a • 2N  (2.30) 


Here  ^ is  the  length  o£  the  one  dimensional  box  and  is  equal  to 
(2N  + 1)D,  where  D is  the  carbon-carbon  bond  length,  which  wc  take  to 
be  1.4  ^ , X is  the  distance  along  the  electron  path  measured  from  one 
end  and  m is  the  mass  of  the  electron. 

Because  of  the  simplicity  of  the  free -electron  orbitals  it  is  possible  to 
Bvibstitute  them  directly  into  the  matrix  components  of  £q.  (2.7)  and  eval- 
uate the  resulting  integrals  directly.  However,  these  integrals  can  be  ex- 
pressed in  a simple  closed  mathematical  form  only  if  the  conjugate  double 
bond  molecules  are  parallel  and  colinear.  For  any  other  orientation,  the 
resulting  answers  can  only  be  expressed  in  numerical  form.  Thus,  it  is 
desirable  to  use  the  London  approximation  which  leads  to  very  simple  re- 
sults with  very  simple  interpretations. 

The  total  energy  of  dispersion  is  the  sum  of  contributions  from  the  sigma- 
sigma,  pi-pi,  and  sigma-pi  electron  interactions; 


E- 


+ p + k 

(KT  FT  ■ C IT 


(2.31) 


The  treatment  of  these  three  types  of  interactions  are  quite  distinct. 


Determination  of  Eg-g- 

A satisfactory  treatment  of  E has  been  given  independently  by 
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London  and  by  de  Boer  in  terms  of  the  polarizabilities  of  the  chemical 
bonds,  in  the  parallel  and  perpendicular  directions,  06  n and  od  ^ 
respectively.  The  mean  bond  polarizability  is  oi  = (06,,  -t  2 ocj  j /3 
and  d ~ ' CC  - '^^')/ od  is  a measure  of  the  electrical  assymetry  of  the  bond. 
In  the  London  or  de  Boer  treatment,  the  energy  of  the  principal  transition 
(h  y or  h depending  whether  the  perturbation  is  parallel  or  perpen- 

dicular to  the  bond)  plays  an  important  role.  We  assume  that  this  energy 
is  equal  to  the  bond  ionization  potenUal,  I . Thus  I = hV^I  = h . 

For  calculational  purposes  we  assume  that  I = 100,  000  cm  ^ = 12.4  e.  v. 
With  these  assumptions,  the  London  or  de  Boer  treatment  gives: 


CTT 


--L' 

^ k 1 


f3 


Ab 


[3 


COS  G - M. 

I< . 


^ K [3 


cos  8 


(2.32; 


S. 


$i"  i Sin  6 cos'^  -J 

•*e  V -V  ' 


^ cos  9 

''a  -^6 


i 

- CCS  0.  - cos  ^3 

k J^a 

B 


J 


5.  F.  London,  J.  Phys.  Chem.  46  , 30  5 (1942);  J,  de  Boer,  Physica  9, 
363  (1942);  also  J.  O.  Hirschfelder,  C.F.  Curtiss,  and  R.B. Bird,  ~ 
“Molecular  Theory  of  Gases  and  Liquids*’,  (John  Wiley,  1954),  p.  969. 
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Here  the  k and  J signify  the  chemical  bonds  in  molecules  A and  B 
respectively  and  the  summation  is  to  be  taken  over  all  of  the  chemical  bonds. 

The  angles  0 ■ and  0 , are  the  angles  which  the  bonds  make  with 

-^6  A ^ 

a radius  vector  joining  the  centers  of  the  two  molecules.  The  JD.  and  (pn 

/3 

are  the  angles  which  the  projection  of  these  bonds  make  with  a plane  passing 
through  the  radius  vector.  Values  of  the  bond  polarizabilities  have  been 

estimated  by  Denbigh^,  For  the  aliphatic  C“C  bond  and  for  the  C-H  bond: 

as"  a 


OL 


cm 


s 


6.  4 

2.90 

(2.33) 

(C-H) 

6;  5 

0.  32 

If  the  two  molecules  are  held  at  a fixed  separation  between  their  cen- 
ters and  rotated  through  all  orientationSj,  we  obtain  for  the  averaged  value 
of  E 


(rq-  j 

E 


ccr 


3 I xn 

^ ^ k ( * ‘ 


(2.34) 


The  values  of  the  E q.  g-  as  calciilated  by  Eq.  (2.12)  are  probably  not 
accurate  to  within  a factor  of  two  but  it  is  difficult  to  see  how  substantial 
improvements  can  be  made  in  the  formulation  without  making  detailed  cal- 
culations using  accurate  wave  functions  for  each  of  the  molecular  species. 
This  short  summary  of  the  treatment  of  the  E j-jj-is  included  here  since  we 
need  to  know  values  of  the  E in  order  to  compare  the  calculated  with 

the  experimental  values  for  the  dispersion  energy. 


6.  K.  G.  Denbigh,  Trans.  Faraday  Soc.  36,  936  (1940) 
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Determination  of  Etttt 

In  seeking  the  simplest  determination  of  £ -pj-  , we  used  the  "free- 
electron**  model  for  the  molecular  orbitals  in  Eq.  (2.4).  Then,  we  used 
the  London  method  of  approximating  the  resulting  integrals.  The  combina- 
tion of  the  two  approximations  lead  to  comparatively  simple  expressions  for 
E ^ ^ which  are  easy  to  interpret.  Our  results  are  to  be  compared  with 
the  more  accurate  calculations  of  Coulson  and  Davies  who  used  good 
molecular  orbitals  and  accurately  evaluated  their  integrals. 


Determination  of  Eq-rr 

The  sigma-pi  electron  interactions  have  not  been  considered  by  previous 
investigators  as  a source  of  inter  molecular  energy.  It  is  therefore  necessary 
to  develop  a new  approach.  Consider  the  interaction  between  the  sigma  elec- 
trons in  molecule  A and  the  pi-electrons  in  molecule  B,  The  sigma  elec- 
trons are  localized  in  chemical  bonds.  For  the  pi-electrons  we  use  the 
London  notion  of  monopoles.  The  matrix  components  of  Eq.  (2.7.)  then 
become  the  sum  of  interaction  integrals  between  the  n^  sigma  electrons 
in  the  p-th  chemical  bond  of  molecule  A with  the  k -th  monopole  in  moleC/  . -le 
B.  The  matrix  components  of  Eq.  (2.7)  then  become: 


A 


(2.35) 


If  the  size  of  the  chemical  bond  is  small  compared  to  the  distance,  R , 

pk 

from  the  center  of  the  p-th  chemical  bond  of  molecule  A to  the  k-th 

B 

monopole  of  molecule  B,  then  I /^\tn  ~ can  be  ex- 


panded in  inverse  powers  of  R 


pk' 


Let  the 


‘fcis  be  the  axis 
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of  the  chemical  bond  and  let  the  center  of  the  bond  be  the  origin.  Then  the 
spherical  coordinates  of  the  electron  ”a"  are  (r  , 0 , ^ Similarly 

St  St  (2. 

the  spherical  coordinates  of  the  monopole  are  ( ^ 

Then^ 


(k)| 


L Z u J ■>  " ' ' p'<' 

' \ p\  / 


- p ; 

4 ' 


(2.36» 


The  leading  terms  of  this  expansion  are: 


K',,< 


f? 


pk 


Cci  a.  cos  0 . Sir.  $2^^ 


'pK  " ~ f'Q.  rpK 

+ Sin©  sirO  , cos  j6  Coo  2^  , 

pk 


(2.37) 


=.  ^ 

K'. 

1-^ 


^yo.Xk 


fL  + 


R 


Since  (a)  is  orthogonal  to  ;/).  (d)  when  we  sxibstitute  this  expansion, 

^ ^ J 

Eq.  (2.  31),  into  Eq.  {Z.Z'),  the  tei’m  1/R  , does  not  contribute.  However, 

pK 

the  next  term  does  not  lead  to  vanishing  integrals  and  it  plays  the  leading  role 
in  the  sigma <^i  dispersion  forces.  Thus  to  a good  approximation. 


7.  "Molecular  Theory  of  Gases  and  Liquids",  Eq.  (12.1-25),  p.  842. 
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ah 


^ e r 


R 


pk 


j <p^  <i>.  Ca.)[X,Xp^  + y 2p,,J  d 


(2. 38) 


Now  the  transition  (i,j)  is  associated  with  a particular  chemical  bond  p 
so  that  the  x component  of  the  transition  dipole  moment  of  the  chemical  bond 
is  given  by 


5y 


(2.39) 


Then  raaking  use  of  the  equivalence  of  the  electrons  and  Eq.  (2.  5)  it  follows 
that 


K d 


The  matrix  components  can  therefore  be  written 


(2.40) 


Rpk 


(2.41) 


Now  we  must  make  use  of  the  special  properties  of  bond  functions  before  sub~ 
stituting  £q.  (2.41)  back  into  Eq.  (2.6).  These  properties  are; 

(1)  Each  transition  i occurs  in  a particular  chemical  bond  so 
that  the  summation  over  p in  Eq.  (2.  41)  is  redundant. 

(2)  If  a particular  transition  i — > j has  a transition  dipole  moment 
in  the  z>direction  parallel  to  the  bond,  it  cannot  have  a transition  dipole 
moment  in  a direction  perpendicular  to  the  bond;  and  vice  versa. 

(3)  The  symbols  i and.  j denote  the  sets  of  quantum  nunnbers  defining 
the  state  of  the  one  electron  bond  orbital.  These  quantum  numbers  are  quite 
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analogous  to  the  set  of  four  fn,  Jt  , m,  cT  ) defining  the  state  of  an  electron 
in  an  atom.  Here  now  the  energy  of  the  state  depends  on  the  quantum  numbers 
corresponding  to  n,  Jl  ^ 'm  | . 

Making  use  of  these  properties  Eq.  (2.6)  becomes; 


1 . 


=.-\ 


'-J  P- 


d b 


pk 


J 


^ ^ \ ✓ 3 ' V 


To  redur-f  this  somm^ilion  still  further,  we  note  that  the  cross-term 


i 


'/'-V  i'-) 

/ > ' If  / I 

Fi’rthermore  the  <len('minator  of  Eq.  (2.40)  depends  only  on  the  absolute 


/ 

V /! 


IS  proportional  to 


jm  - 1 


]. 


value  of  m and  of  m.  and  in  the  summation  both  m.  and  -m.  as  well  as 
1 .1  11 

On  these  accounts  it  follows  that  the  ( ross-terms  vanish 


m and 
J 

and  so 


•m, 

J 


‘-d  "d 


(2.42) 


Now  let  us  assume  that  in  molecule  A there  is  only  on-  principal  transition. 

jn  eH<  h chemical  bond  and  that  the  energy  of  this  transition  is  tj  ‘ = I . 

the  ionization  potential  or  roughly  100.000  cm  ^ or  12.4  e.  v.  Also  the 

9 

parallel  and  perpendicular  bond  polarizabilities  are 


8.  H.  Margenau,  Rev.  Mod.  Phys.  ■ 1 (1939),  see  discussion  on  p.  20. 

9.  See  "Molecular  Theory  of  Gases  and  Liquids",  p.  889,  Eq.  (12.6-l6a). 


2. 17 


(2.43) 


Then,  the  dispersion  forces  between  the  sigma  electrons  of  molecule 
and  the  pi  electrons  of  molecvile  B can  be  written  in  the  form, 


A 


(2.44) 


The  total  sigma-pi  interaction  energy  between  the  two  molecules  is  then 
the  sum  of  E ^ and  E ^ . 


Monopoles  for  Linear  Polyenes 

The  calculation  of  the  monopoles  for  linear  polyenes  is  easily  carried 
out  using  the  FEM  orbitals  of  Eq.  (2.  2')),  assuming  that  the  polyenes  are 
truly  linear  and  that  all  carbon-carbon  bonds  have  the  same  length,  1.4  ^ . 
For  linear  polyenes  the  charge  density  of  the  monopoles  is  given  by  Eq.  (2.8) 

V 

P - ^ € sin  X 

hj  ^ I i (2.45) 

This  has  nodes  »*/h':n.ever  x is  citner  a multiple  of  Jiji  or  a multiple  of 

-£/j.  Thus  the  first  monopole  lies  between  x - 0 and  x =:  .//j  and  each 

of  the  oth.-r  monopoles  lie  between  successive  nodes.  If  y is  the 

k 

X coordinate  rt  the  k-th  node,  then  the  strength  of  the  k-th  monopole  is 


/ 


The  notation 


The  results  for'a  few  polyenes  are  given  in  Table  2,  1, 

(i,  j ) signifies  that  the  monopole s are  for  the  transition  in  which  the  excited 
state  contains  the  j-th  FEM  MC  in  place  of  the  i-th  FEM  MO  in  the  ground 
state . 

In  each  case  the  sum  of  the  monopoles  must  be  zero  because  of  the 
orthogonality  of  the  wave  functions.  The  signs  of  the  monopoles  for  a given 
transition  may  be  changed  simultaneously  without  affecting  the  result,  since 
the  quantity  of  significance  is  the  square  of  the  matrix  element. 

For  a polyene  with  N double  bonds,  the  transition  (N,  N + 1)  represents 
the  largest  contribution  and  is  the  only  one  that  is  listed  for  hexatriene.  We 
shall  speak  of  such  transitions  as  principal  transitions. 


/ 
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•fable  2.  1 Monopoles  for  Linear  Polyenes, 

The  numbers  given  above  the  lines  are  the  strengths  of  the  monopoles 
^^.(k)  given  in  units  of  e.  The  numbers  given  below  the  lines  are  the  dis- 
tances in  % between  the  monopoles.  The  transitions  (ijj.)  are  given  on  the 
left  hand  side  of  this  chart. 

ETHYLENE 

*.C003\  -.C0O3I 


(1.2) 

i7Sa5  — 

ACETYLENE 

(1.2  ) 

flQOO^S.  -/.a«J04<2. 

• « 

f-l.53ii  ■* 

BUTADIENE 

(1.3) 

-sgwf- 

• • • 

(1,4) 

f.n.w  -■4o‘}ts'  -Ji.’in 

0. — ■ •- — — — # 

*-li,959-*-  ♦/.6/93-> 

(2.  3) 

-./o7(.9  +./«7*9  -.4(,TSa, 

(2.4) 

+,3ooio  -.30010  -.30010  +.3ooio 

• 0 ^ ^ # 

♦-a.oi+t -*• 

HEXATRIENE 


- Ol’fCS  *.igo3l  -./toil 

(3,4)  • 
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III.  Results. 

A.  Ethylene 

Ethylene  is  the  simplest  example  of  a molecule  containing  pi-electrons. 
We  shall  now  discuss  the  calculation  of  the  dispersion  energy  for  ethylene. 

We  first  consider  the  energy  arising  from  the  interaction  of  pi-electrons  in 
one  molecule  with  the  pi-electrons  in  the  other  and  shall  compare  the  results 
with  those  obtained  by  Coulson  and  Davies  who  used  the  LCAO  orbitals.  Then 
we  calculate  the  complete  dispersion  energy  by  taking  account  of  the  sigma- 
electrons  as  well  as  the  pi-electrons,  and  we  compare  with  the  experimental 
data. 

The  qualitative  behavior  of  the  ^ ^-fr  interactions  in  ethylene  are 
easy  to  understand  in  terms  of  Eqs.  (2.  6)  and  (2.13).  According  to  Table 
2.  1 there  is  only  one  impcrt9.nt  transition  in  ethylene,  (1,  2)  and  there  are 
only  two  monopoles  associated  with  this  transition.  These  monopoles  being 
equal  in  magnitude  and  opposite  in  sign  correspond  to  a transition  dipole  having 
a dipole  moment 

Aj,' (3.1) 

According  to  Eq.  (2.  30) 

a. 

-II 

G,  =— a(4-i)-  x/0  er^s  = 6-39  e.v.  (3.2) 

For  separations  R between  the  centers  of  the  molecules  large  compared  to 
the  separation  between  the  monopoles,  Eq.  (Z.6)  reduces  to 

^itrr"  ’ hcose* cosset  sinS  cos(si^-  (6J]  ,3  3, 

O _4 

- 3.  Vo  3 ^10  C ergs  ) 


V 


where 


(3.4) 


3.?- 


and  the  angles 


0 


-B 

with  a radial  sector  joining  their  centers  and  and  0 

^ Ij 


are  the  angles  which  the  two  molecules  make 

are  the  azimuthal 

angles  which  these  molecules  make  with  a plane  passing  through  the  radius 
vector. 

At  smaller  separations  Eq.  (3.3)  is  not  quite  correct  because  the 

monopoles  form  a real  dipole  with  a separation  of  1. 7825  % between  the 

charge  elements  as  contrasted  to  an  ideal  dipole  with  txo  separation  between 

the  charge  elements.  The  lack  of  ideality  of  the  dipole -dipole  interaction  can 

be  taken  into  account  by  introducing  arunction  (R,  >9  . » © .q*  a* 

A Ij  a S 

into  the  equation  for  E ^ ^ , 

^TITT  Q^-f-sln  6^Sine^co5(p{^.^^)J  (3.5) 

The  values  of  X (R,  S^'a.  ) are  calculated  using  Eq.  (2.  13) 

A Ij  a a 

to  compute  \ j.Iues  of  different  configurations.  In  the  case  of  ethylene 

where  there  is  only  one  principal  transition,  (1,  Z)  and  only  a single  pair  of 
monopoles  Eq.  (2.  13)  can  be  written  in  the  simple  form: 

G r I 1 I I -j  ^ (3.6) 
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R 


3 • a-a 

are  the  distances  between  the  i-th  and  j-th  monopoles  as 
Fig.  (3.  2);  and  being  the  monopoles  on  one 

Mj  a.nd  the  monopoles  of  the  other  molecule.  Thus  the 


Ra 


R 


J 


where  the  R 
shown  in 
molecule  and 


M.M, 
1 j 


following  results  were  obtained; 

(a)  Parallel  and  Opposed. 


C =:  C 

c = c 


s (%) 


The  results  for  this  case  can  be  summarized  by  the  equation; 

% ey. 

Here  . Y is  a slowly  varying  function  of  the  separation  d {%  ) 
which  should  approach  a constant  value  at  large  separations.  In  Table  3.  1 
the  present  results  (HH)  are  compared  with  those  of  Coulson  and  Davies 
(CD)  who  used  JLCAO  orbitals.  The  values  of  E ^ ^ or  of  V which 
Coulson  and  Davies  calculate  are  uniformly  4/3  of  our  values  (except  for 


the  calculations  at  50  ^ separation). 


3.  3 


Table  3.  1,  Ethylene  Parallel  and  Opposed 


Calculated  Values  of  the  Coefficient 


Y Defined  by  Equation  (3.  V) 


d(^  ) 

^CD 

^ hh/ 

4 

0.760 

1.009 

1. 328 

8 

0.929 

1.245 

1. 340 

10 

0.955 

1. 280 

1.  340 

15 

0.981 

1.  318 

1. 344 

20 

0.988 

1.  315 

1.331 

50 

1.000 

1.079 

1.079 

(b)  Parallel  and  Displaced; 


c=c 


For  this  case,  R 


2 2 

y +d;  cos  coH  3 ~ y ; 

^ ^ R 

0.  Thus  Eq.  (3.5)  can  be  written 


and 


TTV 


01.  i 5 / o ^ 

^ ) If  e.\ 


R"a) 


(3.8) 


Comparing  the  values  of  calculated  from  Eq.  (2.13)  and  the  ex- 
pression for  (3.8)  we  obtain  the  values  of  )(  given  in  Table 

3.2. 


tv. 


/ 


3.5 


Table  3.  2 

^ For  Ethylene  Parallel  Displaced  (Eq.  3.8) 


y (^) 

o 

dr:3A 

d=:8.^ 

0 

0. 6312 

0.9289 

1 

0.4722 

0.9376 

2 

8. 3913 

1. 1527 

3 

0.4918 

1.0052 

4 

0. 8514 

1.0677 

5 

0.9819 

1.3140 

6 

1.0305 

.4625 

7 

1.0451 

. 8475 

8 

1.0471 

.9197 

9 

1.0458 

.9547 

10 

1 . 0420 

.9718 

11 

1.0363 

12 

1.0340 

.9932 

13 

1.0255 

14 

1.0267 

1.0031 

15 

1.0237 

16 

1.0071 

18 

1.0060 

20 

1.0071 

Figure  3.  1 gives  the  values  of  for  d * 3^  and  shows  the  excellent 

agreement  which  is  obtained  with  the  more  exact  calculations  of  Coulson  and 
Davies. 


The  most  characteristic  feature  of  the  curve  for  E^^^is  the  node  which 
occurs  between  y = zX  and  y ~ 3.^  for  the  case  of  d - 3^  . If  the  mono- 
poles  behaved  like  an  ideal  dipole,  then  according  to  Eq.  (3.  7),  the  node  should 
occur  when  y = d/  't  2 or  for  d = 3 .X  , the  nodal  y should  be  2.  12  ^ . 
Actually,  since  the  monopoles  lie  a distance  of  1.78  ^ apart,  the  corrections 

for  lack  of  ideality  are  very  large  for  this  small  separation.  This  is  seen  from 

pC/N-i  U - 3^. 

the  large  value  of  » in  Table  3.  2 for  y ='2A^  Actually  we  find  that  the 
node  should  occur  for  y = 2.43  This  is  in  almost  perfect  agreement  with 

Coulson  and  Davies  who  find  the  node  occurring  for  y = 2.  45  . For  d - 8 

we  obtain  equally  good  agreement  with  Coulson  and  Davies  for  E .p-  as  a 


/ 
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function  of  y.  In  this  case  ws  find  the  node  occurring  at  5.  58  K while 
Cotilson  and  Davies  calculations  lead  to  5.  60  ^ . For  d = 8 ^ the  node 
would  occur  at  y = 5.  66  ^ if  the  monopoles  were  ideal.  The  larger  the 
value  of  d , the  more  closely  does  the  node  occur  at  the  ideal-dipole 
value,  y = . 707  d. 

(c)  One  Molecule  Rotated  In  Plane  of  Other  Molecule  and  the  Radius 
Vector  

c I c 

COS  X (3,9) 

For  R K 7 ^ , the  values  of  V which  we  calculate  are  given  in  Table  3.  3. 
For  this  ease  the  ideal-dipole  approximation  is  very  good.  Coulsan  and  Da- 
vies obtained  the  same  qualitative  behaviour  as  oiir  results  except  that  here 


and  ^ 


B 


According  to  £q.  (3.3) 
for  this  case, 

a I as  y 
^TTfr  R(A)^ 


Table  3.3  X For  Ethylene  Rotated  in  Plane  of 
other  Molecule  and  the  Radius  Vector.  (Eq.  (3.9)) 
Here  R = 7 ^ . 

X Y 

Degrees 

0 1.101 

15  1.080 

30  1 . 042 

45  0.980 

60  1.006 

75  0.882 

90  1.000 
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I again  their  values  are  4/3  of  ours. 


l! 


Complete  Diaper sion  Energy  for  Ethyl^'ne 

The  formulae  for  E jy  ^ and  for  E are  given  by  Eqs. 

(2..32)  and  (2.44)  respectively.  Let  us  perform  the  calculation  for  12 
separate  orientations  of  the  two  molecules  when  their  centers  are  separated 
by  10  . In  the  calculations  of  E y,  and  of  E g- we  use  the  bond  polar- 

izabilities of  (C-C)  ,,  , in  place  of  C - C since  we  do  not  wish  for  the 

aliph 

bond  polarizabilities  to  include  the  contribution  of  the  pi-electrons.  For  the 
C-H  bonds  we  neglect  the  assymetry  of  the  polarizability  tensor  and  use  the 
average  polarizability  in  place  of  oc,,  and  o'.  ^ . This  greatly  sim- 

plifies the  numerical  calculations.  It  is  convenient  to  number  the  monopoles 
(M)  and  centers  of  bonds  (H)  and  (C)  as  shown  in  Fig.  (3.2) 


Fig.  (3.2) 


? 

I 


» 
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For  these  calculations  the  C-H  bond  length  is  taken  to  be  1.07  A and  the 
C=C  distance  is  1.34  X . The  angle  between  adjacent  C-H  bonds  is  100°. 
Let  the  distances  between  the  centers  of  two  C-H  bonds  be  called  R 


between  the  center  of  the  first  C-C  bond  and  the  center  of  C-H 


'H.H.’ 
1 J 


bond,  R_  „ ; between  the  first  monopole  and  the  center  of  a C -H  bond, 

R ; ^ etc.  The  twelve  cases  considered  are  shown  schematically 

M jH. 

in  ^ Figure  3.3.  For  these  cases  the  expressions  for  E fj- q-  is 


(in  e.  V.  when  distances  are  in  ^ 


-3.93  2 I 


(TO- 


>/  j-S 


R 


H.H. 

'■  a 


- o.us  2 rr  [i+i.‘?(3cos  s -i)J 

-s'- S’  V H-  ' 


dr  [^.+  5,9(3005'' 6^  -))] 


- 0.6  af 


R 


c,c, 


2 5 

0,Z  + 0.29  r COS  + COS  e ) 


+ !*#(  (”2  Cos  0^  CobO  -Si’nO^  sio©  )) 


C, 


(3.10) 


Here  ( © . ,0  ) and  ( ©^  . ) are  the  angles  that  the  C = C bonds 

make  with  the  radius  vector  and  a plane  including  the  radius  vector  joining  the 
center  of  the  two  bonds.  Also,  for  the  cases  considered,  the  expression  for 
E is  (ine.v.  when  distances  are  in  ^), 


/ 


3.9 


Fig.  3.  3.  Twelve  configurations  of  two  ethylene  moleciiles.  The 

Bj,  and  represent  four  orientations  of  the  carbon- 

carbon  bonds.  The  other  orientations  differ  from  these  by  rotations 
of  90°  about  the  C-C  axes  as  indicated. 


Molecules  in  same  plane. 


A 


Second  molecule  rotated  perpendi- 
cular to  plane  of  first  molecule. 

Molecules  in  same  plane. 


B 


2 


Second  molecule  rotated  in  plane 
perpendicular  to  first  molecule. 


Molecules  in  same  plane. 


First  molecule  rotated  into  plane 
perpendicular  to  second. 


Second  molecule  rotated  into  plane 
perpendicular  to  first. 


o 

Both  molecules  rotated  by  9C 
into  new  plane. 


First  molecule  in  xz  plane 
second  in  yz  plane. 


“4  -I--—/- 


First  molecule  rotated  by  90 


o 


Second  molecule  rotated  by  90 


o 


Both  molecules  rotated  by  90 


o 
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crir 


r ‘ I 

V. 


( ) 


3 
M 

H.M3 


/a  a a 

-M-3zi^  ' - ’ '3 -";3 

"jN,  IjM.  '5h,M, 


-(or-,\y  iri  C£X4eS  A,  A,^)  5.3 


[( 


.^H.-2c,I 


(2„-2c) 


(3.11) 


Here  the  Z and  2,,  are  the  coordinates  in  the  direction  of  the  radius 

C M 

T«iCtor  from  C,  to  G_:  and  R.,  = 1.783.  For  general 

orientations  of  the  t'wo  molecules,  the  last  term  in  £q.  (3.11)  becomes 
somewhat  more  complicated  but  can  easily  be  worked  out  from  Eq.  (2.44). 

The  results  for  the  12  different  configurations  of  the  ethylene  molecules 
is  given  in  Table  3.4.  In  order  to  obtain  mean  values  of  the  dispersion 
energy  corresponding  to  spatially  averaging  the  molecules  over  all  orienta- 
tions, first  of  all  we  must  average  with  respect  to  rotations  about  the  C-C 
bonds.  This  is  accomplished  by  averaging  the  results  for  and  A^  to 

give  an  average  A;  aversfgie  and  to  get  an  average  B;  average  C^, 

C^,  Cy  and  to  get  an  average  C;  average  Dj^.  Dy  D^*  and  to 
get  an  average  D.  In  the  above  averaging  each  of  the  orientations  is  given 
equal  weight.  The  resulting  four  orientations  of  the  C-C  axes  correspond  to 
the  four  configurations  of  the  H-H  axes  in  the  interaction  of  two  hydrogen 


/ 


1 1 


moleculss  problftm  wh;ch  Evf't*  and  MaTg^nau  srwdj^d.  The  v found  Th-*  .n  • 
case,  the  most  reliable  mean  values  of  »he  .n’e>-ai  *^on  energies  s -iflarned  bv 
suming  the  weights  O.OSS,  0.25,  0.  il5,  and  0.25  (or  the  configurations  A; 

B,  C,  and  D.  Using  these  weights  we  obtained  the  resuJts  rndt;  aied  u Table 
3.4.  Fixst  of  all,  it  is  surprising  how  little  difference  the  orientatjon  of  fbe 


Table  3.4  Energy  of  interartfon  Between 
Two  Ethylene  Molecules  Separated  by  a Distance  of 
10  . /All  energies  are  gi^’en  m un.,t  s of  10  r. . v. ) 


Orientation 

"^trrr 

0 (T 

■■^G-TT 

E, 

to*a] 

^1 

0.91 

1.93 

2.  28 

5.  12 

^2 

0.91 

1.93 

2.  28 

5.  12 

0.20 

0.85 

0.29 

1. 34 

®2 

0. 20 

0.82 

0.28 

1.  30 

0.00 

1.  18 

0.75 

1.93 

^2 

0.00 

1.  18 

0.  75 

1.93 

C3 

0.00 

1. 13 

0.73 

1.86 

C4 

0.00 

1.  13 

0.  73 

1.86 

1^1 

0.00 

0.79 

0.28 

1.07 

0.00 

0.77 

0 . 27 

1.04 

°3 

0.00 

0.77 

0. 27 

1.04 

°4 

0.00 

0.74 

0.27 

1.01 

t 

Mean  Values 

0.  13 

1.04 

0.64 

1.84 

Experimental  Dispersicn 

Energy 

from 

Lennard*' Jones  Potential 

(6,  12) 

4.06 

C-H  bonds  make  in  the  and  E . 

When  one  uses  the  Lennard-Jones  (6,  12)  potential  to  fit  the  experimental 


10.  A.  Evett  and  H.  Margenau,  Phys.  Rev.  90  , 1021  (1953). 
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daita  of  viscosity  of  Ethylene  as  a function  of  temperatiire,  the  force  con- 
stants 6 and  C are  foxind  to  be  0.0177  e.v.  and  4.232  .$  re- 
spectively. If  this  Lennard-Jones  potential  were  correct,  the  energy  of 

6 

dispersion  would  be  equal  to  4 6 C /R  . This  empirical  potential  energy 
curve  therefore  provides  an  experimental  comparison  of  our  average  value 
of  E^  ^ The  fact  that  our  theoretical  value  is  smaller  by  a factor  of 
1.81/4.06  3 0.  45  is  not  surprising  since  most  theoretical  values  of  dis- 
persion energies  are  smaller  than  the  Lennard- Jones  values  by  about  this 
ratio 


t 


11.  See  ^'Molecular  Theory  of  Gases  and  Liquids",  p.  966,  Table  13.3-1. 


/ 
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Acetylene 

The  forces  between  two  acetylene  molecules  are  much  the  same  as 
between  two  ethylenes.  In  acetylene  the  triple  bond  contains  two  pi-bonds  but 
these  pi~electron8  can  be  treated  in  exactly  the  same  manner  as  the  pi- 
electrons  in  ethylene.  The  carbon-carbon  bond  length  in  acetlene  is  1. 207  ^ 
and  the  carbon-hydrogen  bond  length  is  1.060  ^ . Using  the  monopole 
strengths  and  separation  given  in  Table  2.  1 we  obtain  the  skeletal  structure 
of  the  acetylene  molecvile  shown  in  Fig.  3.4.  As  in  the  case  of  ethylene, 
there  is  only  one  principal  transition  which  needs  to  be  considered,  (1,  2),  so 


^ Monopole  Strengths  — ^ 
1.200  e - 1.200  e 


• s s s s 

Hi  Mj  Cj 


-1.134  -0.  768 


0 +0.768  +1.134 


Distances  in 


to  right  of 


Fig.  3.4  , Skeletal  structure  of  acetylene. 
Here  is  the  center  of  the  C-C  bond  and 

Hj  and  are  the  center  a of  the  C-H 

bonds.  The  M^  and  M2  are  the  two 
monopoles  for  the  principal  transition. 


that  the  formulae  for  the  three  different  types  of  electron  interaction  are 
simple.  For  acetylene,  the  transition  dipole  moment  is  ^ debye 

and'  the  transition  energy  is  £ - € . s 8.60  e.  v.  At  large  separations 

E is  given  by  Eq.  (3.3)  where  now  G = 140  e.v.  - .S  . At.«.lXialler 

separations  ^ is  given  hy-Eq.-  (3.^)  where  ..  = 1.537  .X  • The 
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formula  for  E jj-g-  is  given  by  Eq.  (3.  10)  where  now  the  summation  over 
*»i»«  includes  only  1 and  2 and  the  summation  over  "j"  includes  only 
3 and  4.  For  the  purposes  of  computation  for  the  present  paper,  we  con- 
sider only  four  basic  orientations  of  the  two  acetylene  moleciiles,  namely  the 
same  orientations  of  the  C-C  bonds  as  were  considered  in  the  case  of 
ethylene;  A,  B,  C,  D of  Fig.  3.3.  For  acetylene, which  is  a linear  molecule 
it  is  not  necessary  to  consider  the  sub-cases  indicated  by  subscripts  on  the 
A,  B,  C,  and  D.  For  these  cases,  the  sigma-pi  interactions  for  acetylene 
are  given  by 


Here  again  the  Z's  are  coordinates  in  the  direction  of  the  radius  vector 
going  from  to  C^i  and  again  the  formula  for  a general  orientation 

would  have  a last  term  which  is  somewhat  more  complicated. 
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The  results  of  these  calculations  are  given  in  Table  3.  5.  The  mean  values 
are  obtained,  as  in  the  case  of  ethylene,  by  assigning  the  weights:  0.085, 
0.25,  0.415,  and  0.25  to  the  configurations  A,  B,  C,  and  D.  Using  the 

Lennard- Jones  parameters  £•  0.016  e.  v.  and  O’  = 4.221  K 

Table  3.  5.  Energy  of  Interaction  Between 
Two  Acetylene  Molecules  Separated  by  a Distance 
of  10  ^ . (Energies  are  given  in  units  of  10  e.  v.) 


Orientation 

^ total 

A 

1.97 

0.89 

2.  22 

5.08 

B 

0.45 

0.97 

0.  37 

1. 79 

C 

0,00 

1.81 

0.97 

2.  78 

D 

0.00 

0.22 

0.  36 

0.  58 

Mean  Values 

0.  28 

1.  12 

0.  77 

2.  17 

Experimental  Dispersion 
. Energy  from 

Lennard-Jones  (6,12)  Potential  3.62 


determined  from  fitting  experimental  viscosity  of  gas  data,  we  obtain  an 
empirical  check  on  our  theoretical  calculations.  In  this  case  our  theoretical 
dispersion  energy  is  smaller  than  the  empirical  value  by  the  factor 
2.  17/3.62  =•  0.60  which  is  exactly  the  same  as  for  the  case  of  ethylene. 
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B tttadiene 

Now  let  us  consider  the  pi-pi  dispersion  energy  for  Butadiene.  This 
example  is  much  more  complicated  than  either  the  ethylene  or  the  acetylene 
cases.  There  are  four  transitions  which  must  be  considered.  They  are,  in 
the  order  of  their  importance:  (2,  3),  (2,4),  (1,  3),  and  (1,4)  . At  large 

separations,  the  (2,3)  and  the  (1,4)  charge  distributions  could  be  approx- 
imated by  a dipole  but  the  (2,  4)  and  (1,  3)  are  essentially  quadrupole  charge 
distributions.  Naturally,  the  positions  of  the  nodes  and  maxima  in  the  energy 
of  interaction  is  difficult  to  describe  in  terms  of  simple  equations.  Neverthe- 
less, the  '•  long-polyene  approximation”  which  considered  presently  gives  a 
surprisingly  good  representation  at  large  separations.  In  this  approximation 
only  the  principal  transition,  (2,  3),  is  considered.  However,  at  small  separ- 
ations there  is  no  satisfactory  way  of  avoiding  the  detailed  calculations.  Let 
us  use  the  symbol  (i  i*  „ j j')  to  indicate  the  transition  (i  , j)  in  molecule 
A and  the  transition  ( i’  , j'  ) in  the  molecule  B.  The  numbers  preceding 
each  term  indicates  its  degeneracy.  Let  us  consider  two  configurations. 

(a)  Parallel  and  Opposed. 

The  results  for  this  case  are  given  in  Table  3.6. 

From  symmetry  considerations  it  is  easy  to  see  that  the  transitions  (11,  34), 
(21,  33),  (21,44)  and  (22,  34)  do  not  contribute  to  the  energy  of  dispersion. 

It  is  seen  that  the  principal  transition  (22,  33)  is  always  the  largest  contributor, 
followed  by  (21,  34)  and  (22,44).  However,  at  large  separations  only  the 
(22,  33)  is  important. 

(b)  Displaced  Parallel  Configuration. 

Here  the  symbols  y and  d have  the  same  significance 
as  in  the  case  of  ethylene.  The  contributions  of  the  various  transitions  and 
the  values  of  are  given  in  Table  3.7  for  various  values  of  y when 

d = 8 .X  . The  variation  of  with  y is  shown  in  Fig.  3.4,  Because  the 

dispersion  energy  is  the  sum  of  a number  of  terms  each  having  different  nodal 
points,  has  a minimum  value  different  from  zero. 
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It  is  interesting  to  compare  the  results  of  the  present  calculations  using 
the  FEM  orbitals  with  the  calculations  of  Coxilson  and  Davies  using  the  LCAO 
orbitals.  The  first  difference  to  note  is  that  the  degeneracy  of  the  terms  is  not 
always  the  same  for  the  LCAO  and  the  FEM  treatments.  For  example  in  the 
parallel  opposed  configiiration,  in  the  LOAO  treatment  the  (22,44)  is  4-fold  de- 
generate whereas  in  the  FEM  treatment  the  (22,44)  is  singly  degenerate.  How- 
ever, in  the  FEM  treatment  the  singly  degenerate  (11,  33)  and  the  doubly  degen- 
erate (21,  43)  have  very  nearly  the  same  values  for  the  matrix  elements.  There 
is  therefore  very  little  physical  effect  of  this  difference  in  the  degeneracies. 
Rather  xmiformly,  the  CouJson  and  Davies  results  give  values  of  -E.j|,  ^ larger 
than  ours  by  a factor  of  5/3.  Such  a difference  could  arise  through  the  particu- 
lar empirical  value  of  the  resonance  integjral  which  Coulson  and  Davies  used  or 
it  could  result  from  errors  in  the  FEM  approximation  or  in  honest  differences 
between  the  LCAO  and  the  FEM  approximations. 

Hexatriene 

For  hexatriene,  we  consider  only  the  principal  transition,  (33,44) 
whereas  Coulson  and  Davies  cqnslder  a large  number  of  transitions.  Table  3.  8 
compares  o\ir  results  with  those  of  Coulson  and  Davies  for  the  parallel  opposed 
configuration.  Uniformly,  the  Coulson  and  Davies  values  are  5/3  of  ours. 

Table  3.8  -E^p,  (e.v. ) for  Hexatriene 

in  the  Parallel  Opposed  Ckinfiguration. 


^ (^) 

Our  Results 

Coulson  and  Davies 

4 

6.04  X 

10'^ 

11. 35  X 

10-" 

8 

3.  57  X 

10-^ 

6 . 07  X 

10-' 

15 

1. 43  X 

10"^ 

2 . 42  X 

10-^ 

20 

2,  86  X 

10-5 

4 . 85  X 

10-' 

30 

2.75  X 

10-^ 

4 . 70  X 

10-' 

50 

1. 35  X 

10“" 

4.01  X 

10-' 
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Approximate  Treatment  of  Long  Polyenes 

L.et  us  develop  an  approximate  method  for  treating  the  dispersion  ener- 
gy between  linear  polyenes  for  the  case  that  the  number  of  double  bonds,  N,  is 
large  and  the  separation  between  the  two  molecules  is  large  compared  to  the 
length  of  the  moleciile.  Consider  the  energy  arising  from  the  interaction  of 
pi-electrons  with  pi-electrons,  , the  energy  arising  from  the  inter- 

action of  sigma-electrons  with  sigma-electrons,  Eg-j.  , and  the  energy  aris- 
ing from  the  interaction  of  sigma-electrons  with  pi-electrons,  E . 

We  include  only  the  principal  transition  for  the  pi-electrons  so  that 
the  - . transition  charge  density  is  given  by: 


^ e 

A 


. NTrX  . Cw+()T)'X 

Sin  — - — Sin  — : 

X i 


(3.13  ) 


Using  a familiar  trigonometric  identity  this  may  be  rewritten: 


'fz  e 


cos 


CoS 


(3Ni-\)  IT  X 


(3.14  ) 


For  large  N the  second  term  in  the  brackets  is  a rapidly  oscillating  function 
so  we  make  the  approximation  that  it  may  be  neglected.  Our  approximate 
charge  density  is  then  given  by: 


T2  e 

?-  ~r 


cos 


Trx 

1 


(3.15  ) 


In  using  the  London  approximation  we  then  obtain  a "dipole"  and  it  may  easily 

be  seen  to  consist  of  a positive  and  negative  charge  of  magnitude 
6 6 

Q (t)»-Q  (-)  - ('fa /It)  C separated  by  a distance  of  S.  A/ IT  . The 

\n  + I . A N / 2. 

dipole  moment  is  then  (fi'fa  6 H )/  IT  In  Table  3.9  we  compare  the  di- 
pole moment  with  the  dipole  moxxxents  for  the  principal  transitions  for  ethylene, 
butadiene,  and  hfxatriene.. 


r 


Table  3.  9. 

Comparison  of  Principal  Transif  ion 
Dipole  Moments  with  Asymptotic  Val\ie 


Molecule 

Ethylene 
Butadiene 
Hexatriene 
Asymptotic  Value 


N, N+1 

Transition  Dipole  Moment 

O. 25473  ei 
0.27512  e I 
0.28073  e I 
0.28658  e 1 


It  is  seen^^that  for  ethylene  the  deviation  is  not  great  and  that  the  asyunptotic 
value  is  rapidly  approached  which  suggests  that  the  approximation  is  a good 
one.  We  now  make  a further  approximation  that  the  interaction  of  two  such 
dipoles  may  be  given  by  the  interaction  of  two  ideal  dipoles  as  in  Eq.  (3.  3). 
Finally,  making  use  of  the  fact  that  the  length  of  the  molecule  is  (2N  + 1)  x 1. 4 A 
we  obtain  the  following  result  for  the  dispersion  energy  arising  from  pi-pi 


interaction: 


f = -0.l"/003(3Ntl) 
TTTT 


[_;i  cos  G^co^ej^+sin©^  sinQj,  CDS  j 


(3.16) 


where  R,  the  separation,  is  measured  in  K . Averaging  over  angles  introduces 
a factor  of  2/3  and  gives  the  result: 


E ■“  ~0, 0^335 
irtT 


(aN-n)' 


(3.17  ) 


It'  is  seen  that  i>  proportional  to  the  fifth  power  of  the  length  of  the 

polyene,  a result  alsQ  reported  by  Coulson  and  Davies.  However,  their  coeffi- 
cient is  somewhat  larger.  Thus  for  0^=  their  coefficient  is 

0.246  compared  to  the  above  value  of  0.1400  3. 

In  treating  E we  use  Eq.  (2.  32);  treat  all  bonds  as  if  they  were 

at  the  center  of  the  molecule;  and  average  over  all  orientations  of  the  molectiles. 
The  averaging  over  all  orientations  has  the  effect  of  making  all  of  the  terms 
with  ^ *8  vanish  so  that 
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3.2Z 
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(3.  18) 


Here  the  sums  are  to  be  taken  over  all  of  the  chemical  bonds  in  the  two 
molecviles.  Since  a polyene  with  N double  bonds  has  (2N-1)  C-C  bonds 
and  (ZN-fZ)  C-H  bonds,  it  follows  that 

2 

(c.V.  ) (3.19) 


— ^ _ 3.^1  (‘In-*-/) 


or<r 


(?' 


We  have,  of  course,  regarded  all  double  bonds  as  single  bonds  in  this  calcu- 
lation inasmuch  as  we  do  not  wish  to  include  the  pi-electrons  in  this  expression. 
It  is  seen  that  E varies  approximately  as  the  square  of  the  length  of  the 

molecule  so  that  for  sufficiently  long  molecules  it  is  negligible  compared  with 


£ 


iTtr 


In  order  to  obtain  a simple  form  for  the  average  value  of  E g.  ^ for 

long  polyenes  at  large  separations,  let  us  assume:  (1)  that  each  of  the  chem- 

, -24  3 

ical  bonds  has  a mean  polarizability  OC  = 0.645  x 10  cm  and  we  ignore 
the  difference  between  the  perpendicular  and  parallel  polarizabilities  as  well 
as  the  difference  between  the  polarizabilities  of  the  C-C  and  of  the  C-H  bonds. 
(2)  We  assume  that  the  centers  of  each  of  the  chemical  bonds  coincides  with 
the  center  of  the  molecule.  Let  h = distance  between  the  monopoles  and 
and  R be  the  distancesfrom  the  center  of  the  molecule  to  the  plus  and  minus 
monopoles  respectively.  The  geometry  is  then  as  shown  in  Fig.  3. 5. Then: 


Fig.  3.5 
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[[4n  + i j oi 


■<rTr 


4' 


'N-M 


•]L 


i ?.  a 
R_-  h 


(3.20) 


Here  the  factor  2 in  the  numerator  is  due  to  the  fact  that  the  sigma-pi 
energy  is  the  sxim  of  contributions  from  the  sigma  electrons  of  A interact- 
ing with  the  pi's  of  B and  the  sigmas  of  B interacting  with  the  pi's  of  A; 
and  4N  -f  1 is  the  nximber  of  chemical  bonds  in  one  molecule.  From  the 


geometry  of  Fig.3J5  .it  follows  that 

1?^  R^-'-  - Rli  60S  X 

R’  “ + Rh  Goa  X. 


(3.21) 


Substituting  these  expressions  into  Eq.  (3.  20  ) and  expanding  the  resulting 
equation  in  powers  of  h/K  , it  follows  that  the  term  in  1/R^  is  the  leading 
term  and  it  is  given  by 

C3Ca  + h a 

E ' ~ re — ^ ^ ) 

L'^  j J ^ (3.22) 


CTTT 


Now  remembering  that  ) = 19. 17/(2N+1)  e.  v. ; 

I X 12. 4 e.  V.  i h X 2.  8 (2N4’l)alO  ® / TT  ; and  q,  x *(2.  e/  IT 
follows  that 


(Ni-  l an) 


( 1+  3 cos^  X ) 


it 


(3.23) 


Averaging  E over  all  orientations  of  the  two  molecules  corresponds  to 

averaging  dos^  X«  Since  the  spatially  averaged  value  of  C0S^=  , it 

follows  that 


■O-TT 


1.49  (4N-tO  (3lH-h)' 
(n+  / aT)  R ^ 


(3.24) 


3. 2^1 


We  see  tli*t  E varies  approximately  as  the  ctibe  of  the  length  of  the 
molecule  so  that  its  significance  for  large  molecules  is  greater  than  E 
but  less  than  E^^,  . 

, 6 

In  Table3«10'we  present  the  coefficients  of  I/R  for  the  three  dispersion 
energy  terms  for  various  values  of  N , It  is  interesting  to  compare  the  re- 
sults for  N equals  one  with  our  previous  result/,  forethvlene 


6 Table  3.10 

Coefficients  of  I/R  for  Differem  Polyenes  Using 
*<Long  Polyene  Approximation^  The  Coefficients 
are  given  in  Unit*,  of.  e.  v. 


N 

-Err 

-E  R^ 

(TU 

-E  R^ 

irir 

-E  R^ 

total 

1 

98 

89 

23 

210 

2 

317 

512 

292 

1120 

3 

661 

1560 

1570 

3790 

4 

1139 

3520 

5510 

10160 

6 

2440 

8800 

34700 

45900 

8 

4260 

26200 

133000 

163500 

10 

6570 

50400 

381000 

438000 

which  are  given  in  Table 

3.11.  It  may  be  concluded  that  the  long  polyene 

approximation  is  fairly  satisfactory  even  in  the  case 

of  ethylene. 

the  smallest 

polyene . 

Table  3.11 

Comparison  of  **Liong  Polyene  Approximation**  with 

Accurate  FEM  Calculations  for  Ethylei^ 

. The  Co- 

efficients  are  given  in  units  of  b.  v.  -K  , 

, and  have 

the  Vklues  Approached  when  the  Molecules  are 

Separated  by  Large  Distances. 

6 

6 

6 

-Ea-r» 

* ,R 
total 

Accurate 

FEM 

104 

6«i< 

13 

181 

*'LrOng  Polyene  Approx? 

98 

89 

23 

218 

3. 


B»n«ene 

12 

The  FEM  MQ's  for  benaene  which  we  shall  regard  as  having  a cir- 
cular free  electron  path,  are  given  by: 


^0  '■  fair 

’"nc  itr 

~ " J=- 

irr 


Their  energy  is  given  by 

Srr,C^ 


COS  n G 
sin  lO  6 


(325) 


n - I , a,  3 


n = o,  I , a , 3 


(3.26) 


Here  & is  the  angular  coordinate  of  the  electron  and  C ia  the  circumfer- 
ence of  the  ring.  We  assume  that  all  the  carbon-'carbon  separations  are. 

1,4  .X  so  that,  roughly  speaking,  C'*  i < 1 . 4 ^ . Note  that  for  n greater  than 
.vero  the  molecular  orbitals  are  doiibly  degenerate. 

In  the  ground  state  of  benzene,  the  six  TT  electrons  are  located  in  the 
three  orbitals  4 1%  • 

In  treating  benzene  we  shall  consider  only  the  principal  transitions 
which  result  from  transitions  of  'n  electrons  f.rom  either  the  1 c or  the 
1 s orbitals  to  either  the  2 c or  the  2 s orbitals. 

Substituting  the  FEM  MO's  as  given  by  (3.25)  we  have  the  following 
expressions  for  the  transition  charge  densities: 


0.0 

“is  •,  3 C 

sii^e 

C.05  aG 

P 0 

- •‘la;  as  ” 

<3  e 

-;j  SlO  e 

3in  30 

■'g  •’  1C,  i ao 

iri  e 

CO5  0 

cos  i© 

M ^IC-,  35 

i u e. 

cos  e 
T- 

sin  a.® 

(3.2.75- 


13.  H.  Kuhn,  Helv.  Chim.  Acta  (1948),  1780  (1948), 

J.  Chem.  Phys.  1^,  840  (1948),  IT,  1 198  (1949);  Bayliss, 

J.  Chem.  Phys.  287  (1948>;  Simp««n,  J.  Chem.  Phys. 

16,  1124  (1948);,  Ruedenberg  and  Scherr,  J.  Chem.  Phys. 
jr,  1565,  (1953);  Scherr,  J,  Chem.  Phys.  £1,  1582  (1953); 
putt,  J.  Chem.  Phys.  17,  484  (1949),  21,  1597  (195^ 
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The  energy  associated  with  these  transitions  may  be  deduced  from 
Eq,  (3.26)  and  has  the  value: 


G 
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- e. 


7. 

3 h 
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(3.28) 


In  calculating  the  monopoles  it  is  helpful  to  note  that 


- (Ro  tation  by  90°)  x 

§ = (Rotation  by  90°)  x 
a b 


(3.29) 


so  that  we  need  carry  out  calculations  only  for  9^  and  each 

chairge  density  region  we  locate  the  effective  charge  by  calculating  the 
first  moment,  i.e.  , we  obtain  the  mean  values  of  x and  y and  find  that 
in  general  the  monopoles  are  not  located  on  the  ring  but  are  inside  it. 
Thus  we  obtain  the  monopoles  which  are  illustrated  in  Figure  3.6.  The 
position  of  the  atoms  relative  to  the  monopoles  is  of  no  importance  in  the 
FEM  model. 
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2.^  aoT  

a = ( Is  ; 2c} 


-,3col0  +.300/0 


300/0  + 500/0 

<■■-1.5150  — =► 

b = (Is  ; 2s) 


t.Oilif  -.063/3’ 


c = (Ic  ; 2c)  d = (Ic  ; 2s) 


Figure  3.6.  Monopoles  for  Principal  Transitions  for  Benzene 
(Units  are  Angstroms  and  electronic  charge,  e). 
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Let  us  calculate  the  dispersion  energy  between  two  benzene  molecules 
arising  from  pi-pi  interaction  for  the  orientation  in  which  the  two  benzene 
rings  are  parallel  and  facing  each  other  and  are  separated  by  a distance,  d 
(see  Figure  3.  7). 


Figure  3.7 

Parallel  Opposed  Configuration 
for  Benzene. 


Let  us  denote  the  energy  arising  from  the  interaction  of  monopole  dis- 
tribution (la;  2c)  = a with  monopole  distribution  (Is  ; 2s)  = b by  the  symbol 
(a,b,)  and  so  forth.  From  symmetry  considerations  we  see  that  the  terms 
(a,  b),  (a,  c),  (b,  d)  and  (c,  d)  are  zero.  The  terms  which  contribute  are 

(a,  a)  = {c^c),  (b,  b)  - (d,  d)  and  (a,  d)  - (b,  c),  the  last  two  being  each 

two -fold  degenerate.  It  is  then  a straight  forward  matter  to  calculate  the 
different  energy  terms  and  the  results  are  presented  in  Table  3.12-  In 
Table  3.12  the  term  values  are  multiplied  by  their  degeneracy  to  give  their 
contribution  to  the  total  dispersion  energy. 
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Table  3,  12 

- Dispersion  Energy  for  Benzene  for 

Parallel  Opposed  Configxiration  (in  e.  v. ) 


d fj) 

~2(a,  a) 

'•2ib.b) 

-4(a,  d) 

TTir 

5 

3.07  X 10~^ 

1. 17  X lO'^ 

2. 31  X 10“^ 

6.  55  X lo"^ 

0.  70  5 

10 

6.  19  X 10"^ 

?. 36  X 10"^ 

4. 71  X 10"^ 

1. 33  X lO'  ^ 

0.916 

15 

5.  7Z  X 10“^ 

2. 18  X 10'^ 

«/.  Q 

4.  36  X 10 

1. 23  X lO"^ 

0.965 

25 

2.  74  X lO"^ 

1.04  X 10"^ 

2.09  X lO"^ 

5.  87  X lO"^ 

0.987 

50 

4.  34  X 10‘ 

1. 65  X 10‘^ 

3. 30  X lO"^ 

9.  29  X lO"^ 

1.000 

To  a good  approximation,  ^ varies 

inversely  as  the 

sixth  power  of  the 

separati 

on  so  that  it  is 

convenient  to  write  it  in  the  form: 

E 

‘-TTir 

\45  i 

-A''’) 

(3. 30) 

where  ^ is  a slowly  varying  function  with  values  given  in  Table  3.  12. 

13 

Davies  has  made  a rather  extensive  study  of  the  TT-TT  interaction?of 

two  benzene  molecules  using  the  LCAO  method.  His  results  agree  well  with 

the  values  given  in  Table  3,  12  except  that  they  are  larger  in  the  ratio 

228/145  =:  1.  57.  Davies  also  considered  the  variation  of  the  IT-IT  dispersion 

forces  with  angles  and  with  various  types  of  displacements.  He  concludes 

that  the  G’-TT  and  O'- G forces  must  be  much  larger  than  the  TT-TT  forces  in 

order  to  account  for  the  properties  of  crystalline  benzene.  As  a matter  of 
14 

fact  de  Boer  succeeded  in  explaining  the  heat  of  sublimation  of  benzene 

without  ever  considering  TT  electrons. 


13,  P.  L..  Davies,  PH.D.  Dissertation,  King's  College,  University  of 
London  (1949). 

14,  J.  de  Boer,  Trans,  Faraday  Soc,  32,  10  (1936). 
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APPENDIX. 

Polarizability  of  the  Polyenes. 


The  longitudinal  polarizability  of  the  polyenes  arising  from  the 

^ electrons  may  be  calculated  using  equations  (2.44)  and  (2.40).  For  this 

calculation  we  shall  not  assume  the  polyene  to  be  linear.  Rather  let  us  con- 

o 

sider  all  C-C  bonds  to  be  120  and  that  the  polyene  is  in  the  configuration 

o 

(see  Fig.  AJ)»  Because  of  the  angle  of  120  , equation  (2.40)  must  be 
written 


Fig.  A.l 


- /-  A\^/- 


Trans -Polyene  for  Ns3.  The  extra  bonds 
correspond  to  the  free  electron  path  and 
the  dashed  line  represents  the  axis  at  the 
polyene. 


= vr  e (A.l) 

Inserting  the  free -electron  functions  (2.  29)  for  the  transition  ^ ^ I 
for  an  N'polyene  gives  the  result: 


■■  J...  » / 

(,ZhJ 


It  then  follows  from  (2.44)  that  the  polarizability  is: 


o(.  = /•  8l  ^ /o 

N 


-zs 

ZN  i-\ 


C rr\' 


(A.  2) 


(A.  3) 


Results  using  this  equation  are  coznpared  in  Table  A.  1 with  values  obtained 

a 

by  Davies  who  used  LCAO  molecular  orbitals. 

For  large  polyenes  (A. 3)  becomes: 


= /.(?/ X (42/^ +/y  cm 


H 


(A.  4) 


* P.  L.  Davies,  Trans.  Faraday  Soc.  789  (1952) 

/ 
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BO  that  the  polarizability  is  proportional  to  the  cube  of  the  length  of  ihe 
molecule.  Davies  obtains  nearly  the  same  result; 


/«  <■  c ■*  V 


(A . 5', 


Table  A.  1 

Polarizabilities  of  the  N-Polyencs 

N 

H-H 

Davies 

1 

38  ' 10  ^^cm^ 

■JC  ,r-25  3 

35  10  cm 

2 

209 

204 

3 

596 

596 

4 

1290 

5 

2370 

2400 

In  the  case  of  ethylene  it  is  more  significant  to  calculate  the  polar* 
iaability  along  the  axis  of  the  double  bond.  From  Fig.  2 we  see  that  the  di- 
pole moment  integral  (2.40)  must  be  written:  / 


.r  J . _ 


/ 

Fig.  2 Ethylene 


-- 


.)  rf 
4’  ■ ^ t 

« J 

•5.4,/? 


’ y y.  ti  /.  A . 6 

■i  r -A  - 

* y ' ^ r.  ^ ' 

It  is  ihep  found  that  C-'n"  , This  compares  well 

with  the  bond  polarizabilities  of  Denbigh  (Trans.  Far.  Soc.  36,  936  (1940)}^ 

-25  3 

who  found  that  for  a double  bond  is  28.6  10  cm  and  for  a 

-•JR 

The  difference  which  is  presumably  the 


single  bond  is  18.  8 ''  lO'^^cm^ 


-25  3 

contribution  of  the  pi-nlnctrons,  is  9.8  ^ 10  cpi  , which  agrees  well  with 

-25  3 

our  value  of  9.  1'  10  cm  . 


BEST  AVAILABLE  COPY 


/ 


med  Services  I eohnical  Information  figency 


Because  oi’  our  limited  supply,  you  are  requested  to  return  this  copy  WHEN  IT  HAS  SERVED 
YOUR  PURPOSE  30  that  it  may  be  made  available  to  other  requesters.  Your  cooperation 
wilt  be  ai^eciated. 


TiRETSSd  for  any  purpose  other  than  !k  connection  with  a definitely  related 

GOVERNMENT  PROGUHE^tENT  OPERATION,  THE  U.  S.  GOVERNMENT  THEREBY  INCURS 
NO  RESPON^HEILTrY,  NOR  ANY  OBUGATION  WHATSOEVER:  AND  THE  FACT  TEAT  THE 
GOVERi^MEKT  MAY  HAVE  FORMULATED,  FURNISHED,  OR  DI  ANY  WAY  SUPPLIED  THE 
SAID  SJRA WINGS,  SPECIFICATiONS,  OR  OTHER  DATA  IS  NOT  TO  JBE  REGARDED  BY 
mPL^i^ON  OR  OTHERWISE  AS  IN  ANY  MANNER  LICENSING  THE  HOLDER  OR  ANY  OTHER 


PRAV'Tffvf  SPECIFIC  A.TIONS  wvtwwr  oat  a 


iKignON  OR  OTHERWISE  AS  IN  ANY  MANNER  LICENSING  THE  HOLDER  OR  AN^  OTHER 


PERSON  6A  CCKaroRATTON,  OR  CONVEYING  ANY  RIGHTS  OR  PERMISSION  TO  MANUFACTURE 
Ora  OR  SELL  ANY  PATENTED  INVENTION  THAT"MJIY  IN  ANY  WAY  BE  RELATED  THERETO. 


■ /-s?: 


-■  ’ ' ' 


